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Abstract 

, A recent observation of the 511 keV electron-positron annihilation line from the 

q \ Galactic bulge has prompted a debate on the origin of the galactic positrons re- 

sponsible for this emission. Assuming equilibrium between annihilation and injec- 
tion of positrons in the Galaxy, we investigate the possibility that positrons were 
produced by past gamma ray bursts (GRBs). We compare the positron annihilation 
rate inferred by the observed 511 keV line, with the predicted injection rate due 
r S^ ■ to e + e~ pairs leaking out of individual GRBs and those pair produced by GRB 

photons in the dense Molecular clouds at the Galactic center. We conclude that the 
proposed scenario is marginally consistent with observations, and can reproduce the 
observed normalization of the annihilation line only for higher-than-average values 
of the GRB rate in the Galaxy and/or of the Molecular clouds optical depth for 



pair production. 



Introduction — The observations of the 511 keV emission line from the 
Galactic central region with the SPI camera aboard the INTEGRAL satellite 
[1,2,3,4,5,6] have confirmed earlier estimates of the total flux [7,8,9,10,11,12,13,14], 
and provided some crucial information on the morphology of the annihila- 
tion region, suggesting an azimuthally symmetric Galactic bulge component 
with full width at half maximum ~ 8 degrees, and with a 2 a uncertainty 
of ~ 1 degree [5]. The 511 keV line flux in the bulge component amounts to 
(1.07 ± 0.03) 10" 3 photons cm" 2 s" 1 [6], a flux that can be explained in terms 
of electron-positron annihilations via positronium 

T^-^-lO^yr- 1 . (1) 

This annihilation rate provides strong evidence for a large population of Galac- 
tic positrons, which are difficult to account for with standard astrophysical 
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sources. To explain their origin, numerous scenarios have been proposed (see 
e.g. Refs. [5,15,16]), including production by black holes and pulsars [17], mi- 
croquasars [18], radioactive nuclei from past supernovae, novae, red giants 
or Wolf-Rayet stars [16,19], a single recent gamma-ray burst event [20], cos- 
mic ray interactions with the interstellar medium [21] and stellar flares [22]. 
More recently, new scenarios have been proposed invoking annihilation of light 
dark matter [23,24], decay of axinos [25], sterile neutrinos [26], scalars with 
gravitational strength interactions [26], or mirror matter [27], color super- 
conducting dark matter [28], superconducting cosmic strings [29], moduli de- 
cay [30], Q-balls [31], pulsar winds [32], primordial and accreting small- mass 
black holes [33,34]. 

Some of these scenarios are problematic, while others have a wide range of 
uncertainties. For example, the predicted yields and distributions of positrons 
from radionuclei synthesized in supernovae are only marginally compatible 
with observations [35]. Supernovae produce positrons mainly from the decay 
56 Ni -^ 56 Co ^ 56 Fe, although the fraction of positrons that would escape the 
SN ejecta is poorly known. Different estimates suggest that typically 30%- 50% 
of galactic positrons may be explained by SNe la and massive stars (SNII/Ib 
and WR stars). This, combined with the average SN la rate in the Galaxy, 
implies that the contribution of massive stars to the galactic positrons is sub- 
dominant (see Ref. [35] and references therein), although a new estimate of the 
SN la rate and an appropriate choice of the Galactic magnetic field may rec- 
oncile this scenario with observations [36]. In Ref. [16], a new class of SN, SN 
Ic, interpreted as the result of a bipolar Wolf-Rayet explosion, was proposed 
as a source of galactic positrons, from the decays of 56 Co. The supernovae of 
this class are short and bright and may be associated to GRBs [37], but their 
rate is unfortunately unknown. 

It was also suggested that the positrons may be produced from annihilations of 
relatively light dark matter particles, with mass in the 1-100 MeV range [23]. It 
is essential for this scenario that the dark matter particles are light. The usual 
dark matter candidates, with mass in the (0.1-1) TeV range, are expected 
to produce a large number of electrons and positrons through annihilations. 
However, they would also produce high-energy gamma rays. If one requires 
the electron-positron annihilations to occur at a high enough rate to explain 
the 511 keV emission, the associated flux of gamma rays would exceed the flux 
observed by EGRET in the direction of the Galactic Center by several orders of 
magnitude. However, if dark matter is made up of particles with masses below 
the muon mass, the annihilation produces only electron-positron pairs [23]. 

It was shown by Beacom et al. [38] that even if the dark matter particles 
annihilate only into electron-positron pairs, the internal bremsstrahlung emis- 
sion, inevitably associated with this process, would exceed the Galactic diffuse 
gamma-ray emission, unless the dark matter particle is lighter than ~ 20 MeV. 
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A similar constraint applies to other proposed sources of positrons, including 
the decay of exotic particles, such as axinos [25], sterile neutrinos [26], scalars 
with gravitational strength interactions [26], and mirror matter [27]. Recently, 
the upper limit on the injection energy of the positrons in the bulge, thus on 
the mass of MeV Dark Matter, was brought down to « 3MeV through the 
analysis of in-flight annihilations [39]. 

Here we investigate the possibility that galactic positrons were produced by 
gamma-ray bursts (GRB) which took place in our galaxy (see e.g. Refs. [40,41] 
for review of GRB). This scenario is appealing for several reasons. First, GRB 
are known to exist. Second, jets in GRBs are expected to lose energy through 
pair production. Positron production in GRBs has been modeled and was 
suggested as a way to identify the location of a single recent GRB [20], due 
to pair annihilation at the end of the radiative phase. Here we are instead 
interested in those positrons that escape from the fireball and propagate in the 
bulge, as well as in those which are pair produced by GRB photons in the field 
nuclei of molecular clouds. We stress then that our scenario is substantially 
different from those focusing positrons from heavy nuclei decay in hypernova 
ejecta (see Ref. [16]), despite the possible association between these two class 
of objects. 

Propagation of Positrons — Before deriving the rate at which GRBs inject 
positrons in the Galaxy, we first study energy losses and propagation distance 
of positrons in the Galactic magnetic and radiation field. At the energies we are 
interested in, the positrons propagating into the intergalactic medium suffer 
energy losses mainly via ionization, on a time scale [42] 

r — ~ 1Q7 io^T^ (ioS^) yrs ' (2) 

where Nh is the number density of atoms and 7 = E/m e is the initial gamma 
factor of the positrons with energy E. In the Galactic bulge, N H ~ 10 5 m" 3 [42] 
and the resulting stopping time is r; on i z ~ 10 7 — 10 8 yrs, for positrons with 
7 ~ 10 — 10 2 . This time scale is much longer than the typical interval between 
GRBs in our galaxy (see eq. (8) below); therefore, one can treat the injection 
as approximately continuous. 

During the time ri on iz, positrons diffuse in the Galactic magnetic field. Unfortu- 
nately, little is known about the properties of the magnetic field in the bulge. 
One expects this field to have both regular and turbulent components [43]. 
There is compelling evidence of turbulence in the local Galactic magnetic 
field, where the largest scale of the turbulent component is l ce n ~ 50 pc [44]. 
If this turbulence is present in the bulge, then positrons with 7 ~ 10 2 or 
less, have Larmor radii much smaller than the characteristic size of turbulence 
cells. In this regime, we can write the diffusion coefficient for an electron with 
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energy E as (see e.g. Ref. [45] and references therein) 

D(B ).3 Xl 0-(A) 1/3 (_£_)- ,/3 (^ 2/3 cmV1 . ,3) 

This is a phenomenological formula, which only applies when the turbulent 
component is comparable with, or smaller than, the regular component of the 
magnetic field. In the limit of strong turbulence numerical simulations suggest 
that D(E) oc E. Monte-Carlo simulations and analytical approximations for 
diffusive propagation in different regimes have been studied in Ref. [46]. The 
product D(E)T lon \ z provides an estimate of the (square of the) distance traveled 
by the diffusing particle in the case where energy losses are negligible. The 
distance traveled by positrons before being stopped is 

If the regular component is significant, it can help spread the positrons quickly 
over the entire bulge. For example, a straight filament of magnetic field with 
a ~ fiG strength and a negligible turbulent component can transport an elec- 
tron by a distance as large as cr ioniz away from the production site. Magnetic 
field filaments are known to exist in the galactic center, although their global 
configuration throughout the bulge is not known [52]. Therefore, eq. (4) gives 
a lower bound on the size of the positron cloud from a single GRB. 

GRB Rate - We now turn our attention to the estimate of the number of 
positrons injected in the Galaxy. Gamma ray bursts produce photons with 
energies from about 100 keV to 1 MeV over a relatively short time scale, 
1-100 seconds. Based on observations of the afterglows associated with long 
GRBs, one concludes that they have cosmological origin and the energy of 
each burst should be of order 10 52 — 10 53 ergs, assuming isotropy. However, 
there is convincing evidence that the gamma-ray emission is strongly beamed, 
with a characteristic opening angle Q ~ 10~ 3 . In the case of a uniform beam 
(see below for a discussion of structured jets), the energy emitted per GRB 
is [47] 

E n ~ 5 x 10 50 erg. (5) 
Although there is some variation in the observed energy and the opening angle, 
the combined energy estimate in eq. (5) is fairly robust [47,48]. 

Recent observations have provided a strong evidence that GRB is a supernova- 
like event [37]. Supernovae release an energy Estss ~ 10 53 erg, 99% of which 
is carried away by neutrinos, while the remaining 1% can accommodate the 
energetic requirements of a GRB (assuming some reasonable beaming). 

The short time variability on the scale St ~ 10 msec of observed gamma- 
ray signals suggests that the size of the emission region is small, R < cSt 
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3000 km. However, the generation of such a large energy in a small region 
of space would result in massive energy losses due to e + e~ pair production. 
The corresponding optical depth would be as large as r ~ 10 11 , which is, 
of course, unacceptable. To ameliorate this problem, one assumes that the 
emission occurs from a relativistically moving region that has a large Lorentz 
factor, 10 2 ~ 3 [40,41]. In such a model the optical depth can be reduced to 
r < 1, so that photons can escape from the emission region. In addition, 
t < 1 is consistent with the observed non-thermal spectrum of photons. 

The observed local GRB rate is [48,49] 

Observed = 0.5 GpC-V- 1 - (6) 

Taking into account the beaming effect, one infers the actual rate of GRB, 
which is a factor (4tt/Q) ~ 500/q higher. In the simplest jet models fn is a 
parameter of order one [47]. Therefore, for the total observed rate we get 

i2^ = 250/ n Gpc- 3 yr-\ fn ~ 1. (7) 

Assuming a local density of galaxies uq = 0.01 Mpc -3 , we find the GRB rate 
in a galaxy like our own to be 

= 25 A, 10-V" 1 - (8) 

Of course, this estimate assumes that our galaxy is "average" in terms of the 
GRB rate. This assumption may be wrong in two ways. On the one hand, 
the GRB rate may be related to the number of stars, in which case the ac- 
tual rate is higher than our estimate in eq. (8) because our galaxy is bigger 
than average. On the other hand, the rate may be lower if for some reason 
GRBs occur predominantly in smaller galaxies or in galaxies morphologically 
different from the Milky Way. The afterglows are more frequently detected in 
smaller, blue galaxies [50], but, as emphasized in Ref. [50], this may be due 
to a selection effect: one is more likely to detect an optical afterglow from 
an unobscured galaxy than from a galaxy with a large amount of dust. Lc 
Floc'h et al. [50] argue that the optically dark GRBs may originate from dust- 
enshrouded regions of star formation. Eq. (8) implies a galactic GRB output 
in 7-rays 

This estimate was derived under the simplifying assumption of "uniform" 
beaming, i.e. the GRB jet is assumed to be uniform across the entire opening 
angle. Jets could also be "structured", with higher Lorentz factors near the 
rotation axis (see e.g. Ref. [51]). In this case, the energy per GRB and the ratio 
of the observed to true rate, receive opposite corrections, leaving the total 
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energy output of GRBs unchanged [48]. Eq. (9) thus represents a reasonable 
estimate also for jet structures more complicated than a simple uniform beam. 

Positrons injection rate - We write the total energy in electrons and 
positrons as 

Etot,(e+e-) = C(e+e-)-Etot,7 (10) 

where C(e+e-) is the ratio to be determined. A priori, the value of C(e+e-) is n °t 
known and, in a supernova, a "natural" value would be much greater than one. 
However, a successful model of GRB must explain the observed non-thermal 
spectrum of photons. This implies that the photons are emitted from a fireball 
which is sufficiently transparent [40,41]. Note, however, that even for r < 1, 
e + e~ pairs can be produced by different mechanisms with efficiencies up to 

C (e+ e-)~0(l). 

The most efficient mechanism of positron production is pair production of 
GRB photons in the field of nuclei and electrons of the interstellar medium, 
and we denote by Cpair the fraction of the GRB energy transferred to these 
positrons. The environment at the Galactic center is very different from e.g. the 
environment in the spiral arms. Dense molecular clouds populate the so-called 
Central Molecular Zone (CMZ) which is located in the center of the Galaxy 
and which is about 400 pc in diameter [52]. The distribution of molecular gas 
in the CMZ is far from uniform, with clumps of density 10 5 cm -3 , embedded in 
a lower density ~ 10 3 ' 7 cm~ 3 medium. Even higher densities have been argued 
necessary, typically ~ 10 4 (75pc/i? gc ) 18 cm" 3 , in order for the molecular clouds 
at distance R gc from the Galactic center to survive tidal stresses [53] . 

Gamma-rays from GRB explosions crossing the CMZ produce pairs in the field 
of nuclei and electrons of the CMZ gas, with optical depth r CMZ = o~ pp n CMZ l CMZ 
where a pp is the cross section for pair production in the fields of atoms, both 
electron and nuclear fields, while Icmz an d n CMZ are respectively the size 
and the density of the CMZ. For a photon of 10 MeV, we find for Hydrogen 
a^p ~ 0.3 x 10~ 2 barn/atom, and for Helium a^ p e ~ 1.1 x 10~ 2 barn/atom, 
while heavier elements have even larger cross-sections [54]. We take a pp ~ 0.5 x 
10~ 2 barn/atom as a conservative value of the cross-section for the mixture of 
Hydrogen and heavier elements in the CMZ. Adopting a fiducial gas density 
of 10 4 cm -3 (note that the presence of clumps may imply a much larger optical 
depth) and Icmz = 400 pc, we find an optical depth r ~ £ pair ~ 0.03. 

Thanks to the strong, correlated magnetic fields [52] in the central region, 
positrons produced in CMZ would be ejected into the bulge on timescales 
comparable to the crossing time, which is smaller than the (ionization) energy 
loss time even at these high densities. So most of the positrons can escape 
from the CMZ into the bulge. Let us denote ( CMZ the fraction of the GRB 
that either occur inside CMZ or just outside CMZ, with the jet going through 
the CMZ. Based on the stellar density in the bulge, as well as the observations 
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of the afterglows, one can take C CMZ ~ 0.1, also in view of the fact that despite 
its small size, the Galactic bulge contains around 10% of all the Galactic 
molecular Hydrogen and young stars in the Galaxy. 

The spectrum of positrons depends on the spectrum of photons escaping from 
the fireball. Even assuming that the average energy of GRB photons is above 
1 MeV, we estimate the production rate of these positrons to be: 

r t t £tot,(e+ e -) / Cpair \ / Ccmz \ ( 5Me\/_\ . n 49 -1 (u s 

r (e+ e-) )CT ~ Cpair Ccmz ^ ~ ^3 J ) ) ^ ^ " ( U ) 

The injection rate of positrons then appears to be somewhat lower than the 
annihilation rate required by eq. (1). The present scenario thus has to face 
the same astrophysical uncertainties that make the intepretation in terms 
of e.g. type Ic Supernovae, problematic. We stress, however, that although 
the naive estimate of the optical depth adopted above makes the scenario 
marginally compatible with observations, it may underestimate the actual 
optical depth, which is sensitive not only to the average density, but also to 
the actual distribution of gas. Another possible way of increasing rV e + e -) )Cr is 
to assume a higher-than-average local GRB rate, possibly associated with an 
intense starburst in our Galaxy, as suggested by the presence of a large-scale 
bipolar wind at the Galactic center 



Are there other sources of positrons in this scenario? First, some positrons 
can be pair-produced outside the fireball by the outgoing photons interacting 
with the interstellar medium. Observations support this. For example, GRB 
940217 was observed simultaneously by BATSE, COMPTEL, and EGRET in 
the energy range from 0.3 MeV to 100 MeV. The spectrum exhibits a signif- 
icant break at around 1 MeV [56]. The break point, 1 MeV, is right at the 
threshold of pair production, and is likely to be caused by the e + e~ produc- 
tion, as suggested by Winkler et al. [56]. Furthermore, photons can produce 
pairs by interacting with the fireball photons which are back-scattered by 
the medium [57]. However, it has been shown that typically the number of 
positrons produced by this mechanism is rather small [57]. Another source 
is positron leakage directly out of the fireball. The photon pressure inside 
the fireball, at temperatures ~keV, is much greater than the magnetic pres- 
sure [58], and the mean free path and gamma-factor of positrons is comparable 
to those of photons. Finally, some fraction of the positrons can stay inside the 
fireball and escape much later, with a greatly reduced gamma-factor. These 
positrons stay close to the site of the GRB and can be used to identify the 
location of a recent nearby GRB [20]. All these mechanisms do not appear to 
increase significantly the number of positrons predicted by eq. (11), in other 

WOrds C(e+e") ~ Cpair- 

Discussion - Given the uncertainties associated with this and other scenarios, 
it is important to search for signatures that can help distinguishing GRBs from 



7 



other sources of positrons. First, one can test the model by looking for the 
511 keV line emission from nearby galaxies, in particular the nearby Draco and 
Sagittarius dwarf galaxies. The number of stars, the star formation rates [59], 
and, hence, the rate of the GRBs are much lower in these dwarf galaxies than 
in the bulge of the Milky Way. The expected flux from the dwarfs is suppressed 
by many orders of magnitude due to both the low GRB rate and the larger 
distance. 

Boehm et al. [23] have estimated the emission from the nearby galaxies in the 
case of annihilating dark matter, which depends on the integral of the square 
of dark matter density. They concluded that, since the dwarf galaxies are dark 
matter dominated, a large flux could be expected, in contrast with the GRB 
case. The observation of the 511 keV emission from the dwarf galaxies would 
favor the dark matter scenario. In Ref. [60] the search for a 511 keV line from 
the Sagittarius Dwarf Galaxy (SDG) was reported, and no such emission was 
detected. However, the effective observation time of the SDG was only 80 ks, 
and the limits inferred from the data cannot rule out the MeV dark matter 
as the origin of the galactic positrons. In the GRB scenario, the large nearby 
galaxies, such as M31, could be easier to observe than the dwarf galaxies, 
although the flux from M31 would still be suppressed by a factor ~ 10 4 with 
respect to the Galactic bulge due to the larger distance. 

Integral/SPI observations [5] show the ratio of the contributions from the bulge 
(B) and from the disk (D) in the range 3 < (B/D) < 9. This is consistent with 
what predicted from our scenario, since the molecular clouds in the disk are 
less dense, and the production of positron-electron pairs is less efficient. We 
stress that an eventual observation with higher spatial resolution (e.g. with 
the IBIS instrument, also on INTEGRAL), would clarify whether the 511 KeV 
emission is actually diffuse, or due to a limited number of point sources. In the 
latter case, all scenarios involving a diffuse source of positrons (as in the case 
of light Dark Matter), or those involving several sources in a small region at 
the Galactic center, with positrons subsequently spreading in the bulge due to 
the presence of regular magnetic fields (as the one presented here), would be 
essentially ruled out. The intepretation in terms of production in or around 
objects belonging to an old stellar population would instead be favoured. 

If supernovae of a new type, SNe Ic, are associated with GRBs, the positrons 
are produced both through the pair production 77 — > e + e~ during the burst 
and by decaying nuclei [16]. However, the former mechanism dominates over 
the latter in each supernova that also produces a GRB. This is reflected in the 
rates: we need one GRB per million years to explain the same population of 
positrons that could be produced by 200 SNe Ic per million years, as estimated 
in Ref. [16]. If only a fraction of supernovae produce GRBs, the amount of 
pair-produced positrons in the galaxy exceeds that from decaying nuclei as 
long as the rate of GRB per supernova is greater than 5 x 10~ 3 . 
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To summarize, galactic gamma-ray bursts could be responsible for the pop- 
ulation of positrons inferred by the observed 511 keV line from the galactic 
center, but only for higher-than-average values of the local GRB rate or of the 
values of the optical depth for pair production (in the field of molecular clouds 
nuclei) higher than what the estimate presented here would suggest. Although 
such values are not unreasonable, especially in view of the large uncertainties 
associated with both quantities, the scenario appears to suffer of the same 
difficulties of other astrophysical scenarios discussed in the introduction and 
should probably be regarded as only marginally consistent with observations. 

An alternative solution was recently proposed in Ref. [61], where the authors 
suggested that the dynamics of GRB fireballs could be modified in the in- 
nermost regions of the Galaxy, since remnants of GRBs occurring in the hot, 
low-density medium produced by recurrent starbursts in the CMZ, become 
subsonic before they can form a radiative shell. In the formalism introduced 
above, this scenario would translate into a very efficent conversion of the total 
GRB energy output into e + e~ pairs, i.e. C(e+e~) ~ 0(1). Future observations 
should be able to distinguish between this and other proposed models. 

The authors thank J. Beacom and M. Kaplinghat for useful discussions. G. B. 
thanks the TEP group at UCLA for hospitality. S. P.-R. and S. P. thank 
the Theoretical Astrophysics group at Fermilab for hospitality. This work was 
supported in part by the DOE grant DE-FG03-91ER40662 and the NASA 
ATP grants NAG 5-10842 and NAG 5-13399. 
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